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1 DIS kinematics 



The event kinematics of deep inelastic scattering, DIS, are determined by the negative square of 
the four-momentum transfer at the lepton vertex, Q 2 = — q 2 , and the Bjorken scaling variable, 
x = Q 2 /2P ■ q, where P is the four-momentum of the proton. In the quark parton model 
(QPM), the interacting quark from the proton carries the four-momentum xP. The variable y, 
the fractional energy transfer to the proton in the proton rest frame, is related to x and Q 2 by 
y ~ Q 2 /xs, where \fs is the positron-proton centre of mass energy. 

Neutral current (NC) DIS occurs when an uncharged boson (7, Z°) is exchanged between 
the lepton and proton. In QPM there is a 1 + 1 parton configuration which consists of a single 
struck quark and the proton remnant, denoted by "+1". At HERA energies there are signifi- 
cant higher-order quantum chromodynamic (QCD) corrections: to leading order in the strong 
coupling constant, a s , these are QCD-Compton scattering (QCDC), where a gluon is radiated 
by the scattered quark and boson-gluon-fusion (BGF), where the virtual boson and a gluon 
fuse to form a quark- ant iquark pair. Both processes have 2+1 partons in the final state. There 
also exists calculations for the higher, next-to-leading (NLO) processes. 

A natural frame in which to study the dynamics of the hadronic final state in DIS is the 
Breit frame In this frame, the exchanged virtual boson (7*) is completely space-like and 
has a four-momentum q = (0,0,0,— Q = —2xP Breit ) = (E, p x , p y , p z ), where p Bre%t is the 
momentum of the proton in the Breit frame. The particles produced in the interaction can be 
assigned to one of two regions: the current region if their z-momentum in the Breit frame is 
negative, and the target region if their ^-momentum is positive. The main advantage of this 
frame is that it gives a maximal separation of the incoming and outgoing partons in the QPM. 
In this model the maximum momentum a particle can have in the current region is Q/2, while 
in the target region the maximum is ~ Q(l — x)/2x. In the Breit frame, unlike the hadronic 
centre of mass (7*2?) frame, the two regions are asymmetric, particularly at low x, where the 
target region occupies most of the available phase space. 

2 Current Fragmentation Region 

The current region in the ep Breit frame is analogous to a single hemisphere of e + e~ annihilation. 
In e + e~ —>■ qq annihilation the two quarks are produced with equal and opposite momenta, 
± v /s^/2. The fragmentation of these quarks can be compared to that of the quark struck 
from the proton; this quark has an outgoing momentum —Q/2 in the Breit frame. In the 
direction of this struck quark the scaled momentum spectra of the particles, expressed in 
terms of x p = 2p Breit /Q, are expected to have a dependence on Q similar to that observed 
in e + e~ annihilation @, [3|, || at energy ^fs^ = Q. 

2.1 Evolution of \n(l/x p ) Distributions 

Within the framework of the modified leading log approximation (MLLA) there are predictions 
of how the higher order moments of the parton momentum spectra should evolve with the 
energy scale |J |6| . These parton level predictions in practice depend on two free parameters, a 
running strong coupling, governed by a QCD scale A, and an energy cut-off, Qo, below which 
the parton evolution is truncated. In this case A is an effective scale parameter and is not 
to be identified with the standard QCD scale, e.g. Aj^g. In particular, predictions can be 
made at Qo = A yielding the so-called limiting spectrum. The hypothesis of local parton 
hadron duality (LPHD) 0, which relates the observed hadron distributions to the calculated 



parton distributions via a constant of proportionality, is used in conjunction with the parton 
predictions of the MLLA to allow the calculation to be directly compared to data. 




Figure 1: Evolution of the mean, width, skewness and kurtosis of the ln(l/x p ) distribution in 
the current fragmentation region with Q 2 . Data from e + e~ and ep are shown together with the 
MLLA predictions of Dokshitzer et al || (the full line is Qq = A, the dashed Qq = 2A, and the 
dotted Qq = 3A) and the limiting spectrum predictions of Fong and Webber || (dash-dotted 
line where available.) The overlapping points are different x ranges in the same Q 2 range. 
The inner error bars are the statistical error and the outer error bars are the systematic and 
statistical errors added in quadrature. 

The moments of the ln(l/x p ) distributions have been investigated up to the 4th order M; 
the mean (/), width (w), skewness (s) and kurtosis (k) were extracted from each distribution 
by fitting a distorted Gaussian of the following form: 

— - d ° r oc exp (-k - -s6 - -(2 + k)6 2 + -s6 3 + —kd 4 ) , (1) 
a tot dln(l/x p ) y \8 2 4 V ' 6 24 i' KJ 

where 6 = (\n(l/x p ) — l)/w, over a range of ±1.5 units (for Q 2 < 160 GeV 2 ) or ±2 units (for 
Q 2 > 160 GeV 2 ) in \n(l/x p ) around the mean. The equation is motivated by the expression 
used for the MLLA predictions of the spectra p. 

Figure [l] shows the moments of the ln(l/x p ) spectra as a function of Q 2 . It is evident that 
the mean and width increase with increasing Q 2 , whereas the skewness and kurtosis decrease. 
Similar fits performed on e + e~ data show a reasonable agreement with our results, consistent 
with the universality of fragmentation for this distribution at large Q 2 . 



The data are compared to the MLLA predictions of Ref. ||, using a value of A = 175 MeV, 
for different values of Qq. A comparison is also made with the predictions of Ref. |j| for the 
limiting spectrum (Q Q = A). The MLLA predictions of the limiting spectrum in Ref. || 
describe the mean well. However both of the MLLA calculations predict a negative skewness 
which tends towards zero with increasing Q 2 in the case of the limiting spectra. This is contrary 
to the measurements. The qualitative description of the behaviour of the skewness with Q 2 can 
be achieved for a truncated cascade (Q > A), but a consistent description of the mean, width, 
skewness and kurtosis cannot be achieved. 

It can be concluded that the MLLA predictions, assuming LPHD, do not describe the data. 
It should be noted however that a moments analysis has been performed [K|, taking into 
account the limitations of the massless assumptions of the MLLA predictions; this yields good 
agreement between the limiting case of the MLLA and e + e~ data over a large energy range, 
3.0 < x fsZ < 133.0 GeV. 



2.2 Evolution of the x p Distributions 

Scaling violations are predicted in the fragmentation functions, which represent the probability 
for a parton to fragment into a particular hadron carrying a given fraction of the parton's 
energy. Fragmentation functions, like parton densities, cannot be calculated in perturbative 
QCD but can be evolved with the hard-process scale, using the DGLAP evolution [|TTJ equations, 



from a starting distribution at a defined energy scale; this starting distribution can be derived 
from a fit to data. If the fragmentation functions are combined with the cross sections for the 
inclusive production of each parton type in the given physical process, predictions can be made 
for scaling violations, expressed as the Q 2 evolution of the x p spectra of final state hadrons W2 



The NLO calculations (CYCLOPS) |L3| of the scaled momentum distribution exist for DIS. 

The inclusive charged particle distribution, \jo tot da/dx p , in the current fragmentation 
region of the Breit frame is shown in bins of x p and Q 2 in Fig. |2|. The ZEUS || and HI |TJj data 



are in good agreement. The fall-off as Q 2 increases for x p > 0.3 (corresponding to the production 
of more particles with a smaller fractional momentum) is indicative of scaling violations in the 
fragmentation function. The distributions rise with Q 2 for x p < 0.1 and are discussed in more 
detail below. The data are compared to e + e~ data [[T5[] (divided by two to account for the 
production of a qq pair) at Q 2 = s ee . For the higher Q 2 values shown there is a good agreement 
between the measurements in the current region of the Breit frame in DIS and the e + e~ results; 
this again supports the universality of fragmentation. The fall-off observed in the HERA data 
at low x p and low Q 2 is greater than that observed in e + e~ data at SPEAR |Tj|; this can be 
attributed to processes not present in e + e~ (e.g. scattering off a sea quark and/or boson gluon 



fusion (BGF)) which depopulate the current region [17| , 

A kinematic correction has recently been suggested |19| to the NLO calculation [13] of the 
inclusive charged particle distribution which has the form, 1/(1 + (m c s/(Qx p )) 2 ), where m e s is 
an effective mass to account for the massless assumption used in the fragmentation functions. 
It is expected to lie in the range 0.1 GeV < m cS < 1.0 GeV. The x p data are compared to 
the CYCLOPS NLO QCD calculation incorporating this correction in Fig. 0. This calculation 
combines a full next-to-leading order matrix element with the MRSA' parton densities (with 
Aqcd = 230 MeV) and NLO fragmentation functions derived from fits to e + e~ data ||20[| . The 
kinematic correction allows a more legitimate theoretical comparison to lower Q 2 and x p than 



was possible in earlier publications [21j. The bands represent the uncertainty in the predictions 



by taking the extreme cases of m e g = 0.1 GeV and m cS = 1.0 GeV. These uncertainties are 
large at low Q 2 and low x p , becoming smaller as Q 2 and x p increase. Within these theoretical 
uncertainties there is good agreement throughout the selected kinematic range. ZEUS found the 
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Figure 2: The inclusive charged particle distribution, l/cr tot da/dx p , in the current fragmenta- 
tion region of the Breit frame. The inner error bar is the statistical and the outer error bar 
shows the systematic and statistical errors added in quadrature. The open points represent 
data from e + e~ experiments divided by two to account for q and q production (also corrected 
for contributions to the charged multiplicity from Kg and A decays). The low energy MARK II 
data has been offset slightly to the left for the sake of clarity. NLO predictions multiplied 
by the kinematic correction described in the text. The shaded area represents the extreme cases 
0.1 GeV < m c ff < 1.0 GeV. The upper band corresponds to m e g = 0.1 GeV and the lower band 
m cS = 1.0 GeV. 



kinematic correction describes the general trend of the data but it was not possible to achieve a 
good x 2 fit for m eff over the whole x p and Q 2 range. In contrast HI reported at DIS'99 [[TJ[] that 
a good description of the data could be achieved with a value of m eff = 0.6 GeV. The dashed 
line shows the ZEUS NLO calculation from CYCLOPS multiplied by the kinematic correction 
using m eS = 0.6 GeV. There is not a good description of the data, therefore it is concluded the 
distribution also depends strongly on the parameters used to generate the NLO predictions. 
The uncertainties introduced by this kinematic correction restrict to high Q 2 and high x p the 
kinematic range that may be used to extract a s from the observed scaling violations. 



3 Target Fragmentation 



DIS at low x allows a study of fragmentation in the target region following the initial scattering 
off a sea quark (or antiquark). The description based on MLLA [22] is shown schematically in 



Fig. H where the quark box at the top of the gluon ladder represents the scattered sea quark 
plus its antiquark partner. The MLLA predictions are made up of a number of contributions. 
Contribution C, the top leg of the quark box, corresponds to fragmentation of the outgoing 
quark in the current region. Three further contributions (Tl, T2 and T3), which are sources 
of soft gluons, are considered in these analytical calculations to be associated with the target 



region. It is predicted |22| that the contribution Tl( the bottom leg of the quark box) behaves 
in the same way as the current quark C and so should have no x dependence. The contribution 
T2 is due to the colour field between the remnant and the struck quark, and the contribution 
T3 corresponds to the fragmentation of the rungs in the gluon ladder. Both T2 and T3 are 
predicted to have x and Q 2 dependences which differ from Tl. Both the Tl and T2 contributions 
have been calculated and give particles of momenta < Q/2. The collinear gluons T3, on the 
other hand, generally fragment to particles with momentum >jQ/2. For values of the scaled 
momentum x p < 1.0, the region of phase space is analogous to the current region and has 
contributions mainly from Tl and T2. The parton momentum spectra predicted by MLLA, 
over a range of Q 2 and x, are shown in more detail in Ref . || . In the target region these spectra 
are approximately Gaussian for x p < 1; they peak at a value of x p ~ 0.1 — 0.2 in the range 
of x and Q 2 measured by ZEUS f|, falling to a plateau region for 1 < x p < (1 — x)/x (the 
maximum value of x p in the target region). 
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Figure 3: A schematic of DIS scattering at low x within the MLLA framework. Quark C 
represents the struck sea quark in the current fragmentation region. Tl is the other half of the 
quark box which is in the target region. T2 is the t-channel gluon exchange and T3 the rungs 
of the gluon ladder. 

The distributions in ln(l/x p ) are shown for both the target and current regions in Fig. ^. 



The fitted curves shown are two-piece normal distributions |23[] to guide the eye. In contrast to 



the current region, the target region distribution does not fall to zero as ln(l/x p ) tends to zero. 
Although the magnitude of the single particle density at the peak position of the current region 
distribution grows by a factor of about three over the Q 2 range shown, the single particle density 
of the target distribution, at the x p value corresponding to the peak of the current distribution 
(contribution C is equivalent to contribution Tl), depends less strongly on Q 2 and increases 
by only about 30%. In addition the ln(l/x p ) distribution shows no significant dependence on 
x when Q 2 is kept constant. In the target region the peak position of the ln(l/x p ) distribution 
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Figure 4: The corrected ln(l/x p ) distributions for the target and current regions for the 
1995 data. Fitted two-piece normal distributions are shown to guide the eye. The heavy line 
corresponds to the target region, the light line to the current region. The error bars are the 
sum of the statistical and systematic errors in quadrature. 

increases more rapidly with Q 2 than in the current region; this is consistent with the behaviour 
expected from cylindrical phase space. The approximate Gaussian distribution of the MLLA 
predictions peaking at hi(l/x p ) ~ 1.5 — 2.5 is not observed. The results strongly suggest that 
the target distributions are inconsistent with the MLLA predictions when used in conjunction 
with LPHD. 



4 Rapidity Distributions 

There are predictions, based on LPHD, for the rapidity distribution of charged particles in the 
Breit frame [24|. (Rapidity is defined as Y = In ((E + p z ) j(E — p z )) , where E is the energy of 



the particle and p z is the longitudinal component of its momentum.) It is predicted that there 
is a sharp rise in the charged particle density followed by a plateau with a width proportional to 
ln(Q) as one moves from the current to the target fragmentation region. It is also predicted that 
there will be another increase in the particle density and an appearance of a second plateau, 
with the ratio of the two plateaux being 9/4. This ratio reflects the change in colour charge 
from a dominant quark one, in the current region, to the dominant gluonic one in the target 
fragmentation region. 

The rapidity distributions for charged particles are shown in Fig. [5] A flat plateau is 
observed at low transverse momentum, p t . As pt increases , QCD effects gradually evolve the 
flat plateau into an approximate Gaussian, peaking near zero. This illustrates the nature of 
the Breit frame in separating the current and target fragmentation region. The expected step 
between the current and target region of the rapidity spectra is not observed. Also shown in 
Fig. |5| is a comparsion of the data with predictions from LO Monte Carlo models. In general the 
data is well described. ARIADNE |25| agrees well with the data but LEPTO |26[ has problems 



in describing the high p t data and the introduction of soft colour interactions |27j destroys the 





Figure 5: The rapidity distributions of charged tracks in intervals of transverse momentum 
and Q 2 in the Breit frame. The error bars are the sum of the statistical and systematic errors 
added in quadrature. The arrow indicates the position of the origin of the hadronic centre 
of mass system for the (Q) of the data. The histograms show the predictions of LO Monte 
Carlo models, the solid line is the ARIADNE Monte Carlo and the dashed (dotted) line 
the LEPTO ME+PS [2JJ prediction with (without) soft colour interactions. 
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Figure 6: (a) represents the evolution of the coefficient of correlations k with predominant 
variation in Q 2 for corrected data and MC predictions; (b) shows the same quantity where 
predominantly x varies. The corrected values of (Q 2 ) and (x) are indicated for each plot. The 
inner error bars on the data show the statistical uncertainties. The full error bars include the 
systematic uncertainties. 



agreement with the data. 



5 Current- Target Correlations 

The correlation coefficient k: 

k = a c :V i r 1 cov(n c , n t ) cav(n c ,n t ) = (n c n t ) - (n c )(n t ), (2) 

is used to measure the dependence between charged particle production in the current region, 
n c , and production in the target region, n t , where a c and a t are the standard deviations of the 
multiplicity distributions in the current and target regions respectively. For positive correla- 
tions, k is positive whilst for anti-correlations it is negative. At low Q 2 these correlations are 



sensitive to the BGF process which depends on the gluon density of the proton |28[ . 

Figure [] shows the behaviour of the correlation coefficient n as a function of the average 
values of Q 2 and x. Anti-correlations are observed for all values of x and Q 2 ||29|| . The magnitude 



of k decreases with increasing (Q 2 ) from 0.35 to 0.1. According to the analytic results of 
these observed anti-correlations can be due to the 0(a s ) effects (QCDC and BGF). The 0(a s ) 
kinematics in the Breit frame can reduce the particle multiplicity in the current region and 
increase it in the target region. The magnitude of the anti-correlations increases with decreasing 
(x). According to |28| this can be due to an increase of the fraction of events with one or two 



jets produced in the target region. This behaviour is driven by an increase of the BGF rate, 
due to an increase in the gluon density inside the proton. These observations are, qualitatively, 
consistent with the depopulation effects discussed in section [2~^ . 

In addition, Fig. |6| shows a comparison of the data with various LO Monte Carlo models. 
The ARIADNE model agrees well with the data. The LEPTO and HERWIG J50J predictions 
show the same trend as the data but do not reproduce the magnitude of the correlations. 



6 Event Shapes 



The event shape dependence on Q can be due to the logarithmic change of the strong coupling 
constant a s (Q) oc l/lnQ, and/or power corrections (hadronisation effects) which are expected 
to behave like l/Q. Recent theoretical developments suggest that 1/Q corrections are not 
necessarily related to hadronisation, but may instead be a universal soft gluon phenomenon 
associated with the behaviour of the running coupling at small momentum scales These 
non-perturbative corrections are governed by a parameter a . The scale dependence of any event 
shape mean (F) can written as the sum of two terms: one associated with the perturbative 
contribution, (F) pert , and the other related to the power corrections, (F) pow . The perturbative 
contribution to an event shape can be calculated from NLO programs, such as DISENT |32| . 



At HERA a number of infrared-safe event-shape variables have been investigated |BB[ . Their 



definitions are given below, where the sums extend over all hadrons h (being a calorimetric 
cluster in the detector or a parton in the QCD calculations) with four-momentum p h = {E h , p h } 
in the current hemisphere of the Breit frame. The current hemisphere axis n = {0, 0, — 1} 
coincides with the virtual boson direction. 
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Also investigated was the variable i/fj (over the whole of phase space). This variable represents 
the transition value for (2 + 1) — > (1 + 1) jets of the factorizable JADE jet algorithm for a 
particular event. 

A common feature of all the mean event shape values is the fact that they exhibit a decrease 
with rising Q, Fig. [I]. This is due to fact that the energy flow becomes more collimated along the 
event shape axis as Q increases, a phenomenon also observed in e + e~ annihilation experiments. 




Figure 7: The points are the corrected mean values of C and B as a function of Q. The 
inner error bars on the data show the statistical uncertainties. The full error bars include the 
systematic uncertainties. The full line corresponds to a power correction fit according to the 
approach in [3T]. The dashed line is the perturbative (NLO) predction from DISENT ising the 
value of a q found from the full fit. 



A simple ansatz for the power correction would be (_F) pow — T/Q. However the fits using 
T alone are poor and support the more detailed approach outlined in In this approach 
(F) pow is parameterised as follows: 



pow 



- a s (Q) - ^ - + f + l) a*(Q) 



(3) 



where (3q, K are constants dependent on the number of flavours, ap is a calculable coefficient 
dependent on the observable F, \ii is an 'infra-red' matching scale (/i/ = 2 GeV), -M. ~ 1.14 
is a 2-loop correction (known as the Milano factor) and a® is an universal, non-perturbative 
effective strong coupling below /z/. 

The results of the fit are shown in Fig. || The parameter ao is observed to be ~ 0.5 for 
all event shapes (except the jet rate parameter yjj) , consistent with theoretical expectation. 
However there is a large spread in the values of a s . The theoretically calculated parameter for 
the power corrections for i/fj was a yfJ = 1. This is contrary to the observed need for small 
negative hadronisation corrections [Q. A reasonable fit for i/fj was achieved by using a value 
of a yjJ = —0.25. The extended analysis of the mean event shapes in DIS are consistent with 
the application of power corrections according to |31| though there is still need for further 
understanding. 



> 
o 

« 0.6 

o 

IS 

0.5 
0.4 

0.1 0.11 0.12 0.13 0.14 
a s (M z ) 

Figure 8: Results of the fit to ao,a s for the means of r,B,rc,p and C. The ellipses illustrate 
the la and 2a contours including both statistical and systematic uncertainties. 

7 Summary 

To understand the underlying QCD processes in DIS it is necessary to study the hadronic final 
state. At the current level of understanding, QCD works well and describes the HERA data. 
As the precision of the HERA data improves and the NLO QCD calculations become available 
the framework of QCD is being tested more thoroughly. 
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